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A Josephson supercurrent between two singlet superconducting electrodes separated 
by an extended region of half-metallic Cr02 has been recently observed. Since in a 
half metal only electron pairs with equal spins can carry a supercurrent, there is a 
conversion from singlet to triplet currents at the interface. The conversion mechanism 
remains, however, unclear. Another unresolved issue is the symmetry of supercon- 
ducting correlations in ferromagnets depending on the degree of disorder. In our work 
we address both problems. We introduce a conversion mechanism based on a broken 
spin-rotation symmetry around the magnetization axis of the half metal. We show that 
in the diffusive case the supercurrent is dominated by the product of odd-frequency 
s-wave and even-frequency p-wave pairing amplitudes, whereas in the ballistic regime 
the admixture of other components becomes relevant. We study the entire crossover 
from the ballistic to the diffusive regime. 
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Half-metallic ferromagnets have great potential in the 
field of spintronics as sources of spin-polarized electric 
currents. Remarkably, they show conducting or insulat- 
ing behavior depending on the direction of the electron 
spin. Since only a few such half metals are known, among 
them La2/3Cai/3Mn03 [l| and &O2 [21], their character- 
ization has attracted great attention. Half metals, when 
contacted with other materials such as superconductors, 
can also be used as well controlled test-laboratories to 
study the interplay between different types of orders. 

Recently, Keizer et.al Q reported a Josephson super- 
current between two singlet superconducting electrodes 
(NbTiN) separated by a wide region of Cr0 2 . In a half 
metal only conduction electrons with equal spin can be 
paired, since the other spin species is insulating. Cur- 
rently, the mechanism involved in the conversion process 
between singlet Cooper pairs (||J,) — ||i))/v2 and equal 
spin pairs |tt) at the interfaces between the materials 
remains highly controversial. In the same experiment 
Q a second, seemingly unrelated, observation was made, 
namely the presence of a homogeneous biaxial magnetic 
asymmetry in the Cr02 film (see Fig. [IJ. 

The symmetries of the relevant pairing correlations 
mediating the triplet supercurrent and their dependence 
on the amount of disorder is currently debated. It is 
claimed [!, 0] that the main source of the triplet Joseph- 
son current in diffusive ferromagnets is odd-frequency s- 
wave pairing amplitudes. On the other hand, it has been 
shown that in clean half metals also p-wave triplet pairing 
amplitudes are important [B[ . The samples in experiment 
Q are moderately disordered, and a detailed theoretical 
study of the crossover from the ballistic to the diffusive 
regimes is required. 

Here we report an extensive analytical and numerical 
investigation of the Josephson effect in the hybrid struc- 
ture shown in Fig. [1] and present results for an arbitrary 
concentration of impurities in the materials, thereby clar- 



ifying the symmetry properties of the superconducting 
correlation functions involved. We link the mechanism 
for the current conversion to a broken spin-rotation sym- 
metry around the magnetization axis M in the half metal. 
We argue that the second observation in experiment Q 
of a biaxial magnetic asymmetry can be directly related 
to the conversion mechanism of the singlet supercurrent 
to the triplet supercurrent. 

The mechanism we propose leads to a natural ex- 
planation of several findings of the experiment Q: (i) 
hysteretic shifts of the equilibrium phase difference over 
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FIG. 1: Conversion between singlet and triplet super- 
currents. The setup consists of two singlet superconductor 
banks separated by a half-metallic ferromagnetic Cr02 layer 
with a magnetization vector M aligned along the easy axis 
ei. Spin-rotation symmetry around M is broken by the si- 
multaneous presence of the second easy axis e2. As a conse- 
quence, there is a conversion from a supercurrent of singlet 
Cooper pairs (jsingiet , blue lines) in the superconductors to a 
supercurrent of triplet Cooper pairs in the half metal (jtripiet, 
yellow line). The conversion takes place in the superconduc- 
tors in a layer extending about a superconducting coherence 
length from the interface, as illustrated by the shading from 
blue to yellow. 



the junction depending on the magnetic pre-history; (ii) 
Josephson junctions involving half metals are 7r-junctions 
after subtraction of the hysteretic shifts; (iii) sample-to- 
sample fluctuations in the magnitude of the critical cur- 
rent. 

We show that in moderately disordered half metals the 
supercurrent is carried predominantly by odd-frequency 
s-wave and d-wave amplitudes, multiplied with even fre- 
quency p-wave and /-wave amplitudes. In the diffusive 
limit the supercurrent is dominated by the product of 
the s-wave and the p-wave amplitudes. We find that a 
peak in the temperature dependence of the critical cur- 
rent, previously predicted for clean half metals is a 
robust feature also for disordered half metals. 

The Josephson junction we study, shown in Fig.[IJ con- 
sists of a half metal extending from x\ to x%, sandwiched 
between two singlet superconductors. When a phase dif- 
ference \2 — Xi exists between the superconducting order 
parameters, an exotic form of Josephson effect occurs: a 
singlet supercurrent, jsingiet (blue in Fig.[lJ, is converted 
to an equal-spin triplet supercurrent, j triplet (yellow in 
Fig. [U , within an interface layer extending about a su- 
perconducting coherence length into the electrodes. The 
equal-spin triplet supercurrent flows through the half- 
metallic material, while the singlet part is completely 
blocked. The sum of the singlet and triplet currents is 
constant, obeying the continuity equation. 

The conversion process between the singlet and equal- 
spin triplet supercurrents is triggered by two important 
phenomena taking place at the interface: (i) spin mixing 
and (ii) breaking of spin-rotation symmetry with respect 
to the magnetization axis M in the half metal. Spin- 
mixing is the result of different scattering phase shifts 
that electrons with opposite spin acquire when scattered 
(reflected or transmitted) from an interface Q . It results 
from either a spin-polarization of the interface potential, 
or differences in the wavevector mismatches for spin up 
and spin down particles at either side of the interface, or 
both. It is a robust and ubiquitous feature for interfaces 
involving strongly spin-polarized ferromagnets. 

Broken spin-rotation symmetry leads to spin-flip pro- 
cesses at the interfaces. Its origin is more subtle and de- 
serves special attention. Here we propose a mechanism 
that applies to biaxial ferromagnets with local magnetic 
moments coupled to itinerant electrons. Indeed, the fully 
spin-polarized CrC>2 films used in the experiment Q have 
been shown to be homogeneously biaxial [l(|. At the 
interfaces with the superconductors, the local moments 
are expected to show a certain degree of disorder for two 
reasons. First, the coordination number is lowered by 
the broken translation symmetry at the interface which 
leads to enhanced fluctuations in the direction of the lo- 
cal moments from that in the ordered phase. Secondly, 
structural disorder leads to vacancies in the lattice of mo- 
ments. The exact microscopic distribution of local mo- 
ments at the interface is not important for superconduct- 
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FIG. 2: Breaking of spin-rotation symmetry at the 
interface, (a) Misaligned local moments at the interface in- 
volving a magnetic material with biaxial anisotropy prefer the 
second easy axis C2, leading to an explicitly broken spin rota- 
tion symmetry around the magnetization axis (here parallel 
to the easy axis e±. (b) If the magnetization is pointing ex- 
actly between the two easy axes, spin-rotation invariance is 
possible, (c) However, hysteresis effects usually prevent such 
a situation for any orientation of the magnetization vector. 

ing phenomena, since Cooper pairs are of the size of the 
coherence length £ which is much larger than the atomic 
scale. It is, however, important for the effective interface 
scattering matrix, as it can lead to spin-flip terms if the 
directional distribution of the local moments breaks the 
spin-rotation symmetry around M. 

We propose that the presence of two easy axes with 
relative angle a leads to a non-trivial distribution of the 
moments at the interface, where a majority of the mo- 
ments are pointing along the magnetization axis of the 
bulk half metal but with a significant number of moments 
misaligned but preferably along the second easy axes. For 
example, the distribution in Fig. [2ja) leads to a model 
where the spin-rotation symmetry around M (satisfied 
in the bulk material) is broken near the interface within 
a few atomic layers. 

To quantify the above discussion, we employ a simple 
model that is formulated in terms of an interface scatter- 
ing matrix, which connects incoming to outgoing waves 
in the asymptotic regions [5|, 
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Here, iff and are transmission amplitudes from the 
two superconducting spin bands to the conducting half- 
metallic spin-t band, and #, tf-fj, and are spin- mixing 
angles. This scattering matrix is appropriate for the tun- 
neling limit, in which case its reflection amplitudes have 
unity modulus. The extra phases ±$/4 ensure unitarity 
of S up to linear order in the transmission amplitudes. 
For definiteness, we have chosen the spin quantization 
axis in the half metal along M, and in the superconduc- 
tor such that the reflection amplitude of the scattering 
matrix is diagonal in the tunneling limit. In a system 
with spin rotation invariance around M these two quan- 
tization axes coincide. This is, however, in general not 
the case if spin rotation symmetry around M is broken. 
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The presence of the spin- flip terms tjf e ±I1?iT in the scat- 
tering matrix Eq. ([1]) is a direct consequence of the broken 
spin-rotation symmetry around M at the interface. 

In the following we calculate the Josephson current 
through the junction in leading order in t^, t^, and 
1). As corrections to the singlet order parameter A arise 
in second order in we can neglect the suppression of 
A near the interfaces. In this case Anderson's theorem 
11, JjJ holds and A is also insensitive to impurity scat- 
tering. For simplicity we consider the case of equal gap 
magnitudes in the two superconductors, Aj — |A|e IXj , 
for superconductors j — 1 and j = 2, see Fig. [TJ 

Due to spin mixing at the interfaces a spin triplet (S = 
1, to = 0) amplitude ft j( x ) is developed that extends 
from the interfaces about a coherence length into each 
superconductor, 



f toj (x)=i*\A\e*x>0 j 



(2) 



where £2 n — \/ e\ + |A| 2 . We have separated the influ- 
ence of the interfaces from that of the impurity disorder 
in the bulk materials by introducing the real functions 
(Pq'?(x). The superscript denotes symmetric (s) and an- 
tisymmetric (a) components with respect to /x = cos(# p ), 
where 9 P is the angle between the Fermi velocity and the 
x-axis. In the clean limit, (p^(x) = -5eM e -l*-*il/«slfl 
and^O) = -tp^ j (x)sgn(ij,e n ), whereas = vs/2Cl n , and 
vs is the Fermi velocity in the superconductor. For an ar- 
bitrary impurity concentration the ^-functions are mod- 
ified and must be calculated numerically for each given 
value of mean free path, as described in the Methods. 

Whereas the interface value of / f ° 3 does not change 
with varying mean free path, the interface value of / t s OJ 
increases with decreasing mean free path in the supercon- 
ductor, £s: as until it reaches values comparable 
with the singlet amplitude. Their decay length in the su- 
perconductors decreases, and changes from (^g 1 
in the ballistic limit to \/$s^s/3 in the diffusive limit 
(see Supplementary Fig. 1). For defmiteness, in the re- 
maining discussion we use for the mean free path in the 
superconductors £s = 0.1£o with £o = Vs/2icT c . 

The singlet-triplet mixing in the surface layers of the 
superconductors described above, together with the pres- 
ence of spin-flip tunneling amplitudes, leads to an equal- 
spin (m = 1) pairing amplitude f\\(x) in the half metal. 
We would like to emphasize that it is only the m = 
triplet component derived above that is coupled via the 
spin-active boundary condition to the m — 1 pairing am- 
plitude in the half metal. The singlet components in the 
superconductor, being invariant under rotations around 
any quantization axis, is not directly involved in the cre- 
ation of the triplet in the half metal. A picture of an 
indirect Josephson effect emerges, therefore, that is me- 
diated by the appearance of the m — triplet amplitudes 
in the superconductor. 



In the tunneling limit it is convenient to split the pair- 
ing amplitude in the half metal into contributions in- 
duced at the left and right interfaces: = f^i + /ff2> 
with momentum-symmetric and antisymmetric compo- 
nents 



ftf s (x) = 2ni(t nj t m ^)\A\e^ 



s.a / \ 



(3) 



where Xj — Xj ~ $Ttj — $Uj- Here, like above, we have 
separated the contributions from the interface scattering 
matrix and the contributions from the disorder in the 
half metal by introducing the (real) functions ff a ■ 

The Josephson current results from Eq. (fT6|) in Meth- 
ods, 



Jx = -J c sin(x2 - Xi), 
where the critical current density is given by 



Jc — Jo 
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Here, the current unit is Jo = 47rew^fiV^fT c , Nu is the 
density of states at the Fermi level in the half metal, e 
is the electron charge, Aj — (^TTj^lTj^i)' an( i (" ' ) = 
J dfi - ■■ . The effective phase on each side j is 



Xj = x 3 - (^ttj 



(6) 



Eqs. (ID)-© describe an exotic Josephson effect in several 
respects. Eq. (|6]) is related to the phase dependence of 
the Josephson effect and can be tested for example by 
studying the Fraunhofcr pattern in a magnetic field. For 
a half metal, there can be extra phases that lead to shifts 
of the Fraunhofer pattern lla l. Within our model 
there are contributions #j.f2 — and 2 — $tti ^° ^ ne 
phases that depend on the microscopic structure of the 
disordered magnetic moments at the two interfaces. The 
microstructure can be affected for example by applying a 
magnetic field that leads to hysteretic shifts of the equi- 
librium positions depending on the magnetic pre-history. 
When subtracting the shifts, the junction shows the typi- 
cal characteristics of a 7r-junction, [3] as revealed by the 
minus-sign in Eq. |4]). The possibility to manipulate the 
shifts with an external field yields a way to measure the 
combination + -d^ of the spin-mixing angles. This 
is important as there have been no experiments to date 
that give information about the spin-mixing angles. Fi- 
nally, the critical Josephson current is proportional to 
the spin-mixing angles 1? and to the spin-flip rates of the 
two interfaces. This points to a strong sensitivity of the 
critical Josephson current to interface properties and is 
expected to lead to strong sample-to-sample variations. 
All these findings are in agreement with the experiment 
0- 

We now proceed with a detailed description of the role 
of disorder in the materials. Anisotropic superconducting 
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FIG. 3: Critical Josephson current, (a) In the crossover 
from the ballistic to the diffusive limits the critical current 
J c is monotonically suppressed, (b) In the crossover region, 
contributions to the current from higher partial waves (I > 2, 
d-wave, /-wave, etc.) of the superconducting Green's function 
are suppressed. In the diffusive limit the current (black line) 
is given by a product of s-wave (I = 0) and p-wave (I = 1) 
components (the blue line), (c) For given mean free path, J c is 
exponentially suppressed for junction lengths large compared 
with the coherence length (£ c = vh/2tT in the ballistic limit 
and = VscS /3 in the diffusive limit). J c is given in units 
of JoA\A2/<in. The length unit is £o = (,c(T c ). We assumed 
an anisotropy of both tff and tj,f proportional to \fi\. 



correlations are sensitive to impurity scattering. From 
studies of unconventional superconductivity it is known 
that superconductivity disappears at a critical impurity 
concentration [13j ]. This is however not the case for the 
proximity induced pairing amplitudes studied here. In 
Fig. OH we show results for the critical Josephson current 
as function of the elastic mean free path, normalized to 
Co = vh/^kTc. As shown in Fig. 02a), the critical cur- 
rent is monotonously suppressed for decreasing mean free 
path, from ballistic (left part of the abscissa in the figure) 
to diffusive (right part) transport. The suppression is ex- 
ponential in the diffusive region, with a crossover taking 
place at a mean free path £jj comparable with the clean 
limit coherence length £ c = vh /2ttT. 

In order to better understand the nature of the 
crossover we show in Fig. |3]Jb) the decomposition of the 
Josephson current into symmetry components. We note 
that the current is carried by the product of neighbor- 
ing symmetry components of the functions Tp^ = Ajtpj 
in Eq. (03), i.e. J c = Y^iZo ^c;l,l+x- This can be seen by 
expanding the pairing amplitudes in Legendre polynomi- 
als, fpAn) = YmLq Plit^Pj h leading to terms of the form 
<?2,z^u+i - ^1,2^2,2+1 in E q- where I = 0, 1, 2, 3, . . . 
denotes the s-, p-, d-, /-wave etc. pairing components. 
As all pairing amplitudes involved are spin triplet, the 
even I amplitudes are odd in frequency while the odd I 




FIG. 4: Triplet correlation functions in the half metal. 

The contribution to the Josephson current from the s- 
wave triplet equal spin correlation function always enters 
as a product with the p-wave, with one of the two orig- 
inating from the left superconductor and the other from 
the right. For clean half metals, shown in (a), the p- 
wave component is larger than the s-wave. For more dirty 
structures, shown in (b), the p-wave component is sup- 
pressed compared with the s-wave and the Josephson cur- 
rent is suppressed accordingly. The symmetry components 
are defined as F?(x) = T £ £n>0 </ m (z)}^, Ff(x) = 

l T E £n >o(A t /nj( a: ))st 1 i Plotted in units of iAj \ A\e^ . Am- 
plitudes for fixed frequency give a similar picture. 

amplitudes are even in frequency [l6j]. The sum of the 
first three contributions (red dashed line) , composed from 
the s ■ p (blue), p ■ d (green), and d ■ f (purple) compo- 
nents, amount already to almost the entire current (black 
line). In the diffusive limit, the current is carried almost 
exclusively by the product of the even-frequency p-wave 
and the odd-frequency s-wave pairing amplitudes. The 
crossover region to ballistic transport is characterized by 
an onset of contributions from higher order partial waves 
I > 2. It is clear from the figure, that for £jj > £o the 
diffusive Usadel approximation breaks down. 

In Fig. 03 (c) we show for several mean free paths 
the dependence of the critical current on the junction 
length L. A rapid exponential suppression of the effect 
with junction length is observed in the diffusive region, 
whereas in the moderate disordered region a considerable 
effect is expected for junction lengths up to 5-10 coher- 
ence lengths. 

In Fig. OH we present an analysis of the spatial depen- 
dences of the odd-frequency s-wave and even-frequency 
p-wave pairing amplitudes in the half metal. We show for 
the ballistic case (£h = 10£o) and for the diffusive case 
(£h = Co/10) the pairing amplitudes induced from the 
left and right interfaces. By multiplying the two black 
curves with each other and the two red curves with each 
other, and summing the two contributions, we obtain a 
quantity related to the s • p contribution to the Joseph- 
son current [see Eq. ©]. For ballistic systems the p- 
wave amplitudes are larger than the s-wave amplitudes 
near the interfaces, while the opposite holds for diffu- 
sive systems. Note, however, that the amplitudes are 
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FIG. 5: Non-monotonic temperature dependence, (a) 

The critical Josephson current has a maximum at a low tem- 
perature that for a specific junction length [here L = £o = 
fc(Tc)] depends on the mean free path in the half metal. For 
is > Co the peak position roughly coincides with that of the 
tu = Co curve. The current has been normalized to the zero- 
temperature value. In (b) we show the normalized current 
as function of junction length L for a fixed mean free path 
£h ~ f -5£o- From (a) and (b) it follows that for diffusive sys- 
tems the junction becomes effectively long compared with the 
diffusive limit coherence length, L = Co 3> £d(T); the current 
is then dramatically suppressed [see Fig. |3ja)] and the peak 
is shifted to lower temperatures. 

tied to each other through the following general relation 
between the momentum-antisymmetric and the momen- 
tum symmetric parts: — sgn(e n )iiS,Hdxf^p where 
C^ 1 = l^ 1 + 2\En\/vH ■ In the diffusive limit, there is an 
additional relation, /|^r wave = — sgn(£„)^^9 x /* 1 T wavc . ft 
follows (see Fig.fJ} that the magnitudes of the amplitudes 
differ (their ratio depends on the amount of disorder) 
while the decay lengths of the two are always identical, 
crossing over from the ballistic coherence length £ c = 
vh/ZttT to the diffusive coherence length £<j = 

Finally, we discuss the influence of disorder on the tem- 
perature dependence of the critical current, see Fig. [5] 
We have normalized all J C (T) curves to their zero- 
temperature value. There is a characteristic peak ap- 
pearing at a temperature below ~ T c /2 as predicted for 
ballistic systems in Ref. The origin of the peak is 
the extra factor s\l$? n in Eq. ([5|), that results from the 
odd- frequency pairing amplitudes on the superconducting 
sides of the interfaces being the sources of the equal spin 
correlations in the half metal. 

In Fig.[5£a) we discuss J C (T) for mean free paths in the 
half metal ranging from ballistic to diffusive. Clearly, 
the peak is observed also in the extreme diffusive case 
of £h = £o/100, but is for a particular junction length 
L shifted to lower temperatures with decreasing mean 
free path. Since the peak survives the crossover from 
the ballistic to the diffusive limits, it serves as a true 
fingerprint of equal spin triplet correlations in the half 
metal. The peak is readily measurable experimentally 
when L is of the order of the coherence length at T c . For 
increasing L the peak is shifted to lower temperatures, 
as shown in Fig. EKb) for £h = l-5£o- For long junctions 



the current has a characteristic exponential temperature 
dependence above the peak. 

Turning to Ref. 0, the low-temperature resistivity of 
the Cr02 material of 8.9 fifl cm implies a mean free path 
of Ih ~ 40 nm. From band structure calculations, the 
Fermi velocity for CrC-2 is vh ~ 2.2 x 10 5 m/s [13], imply- 
ing £ ss 27 nm, or £ /l H w 2/3 (for T = 0.3T C this gives 
£ c w 90 nm and £ d » 35 nm). The Cr0 2 of Ref. || is 
therefore moderately disordered, beeing in the crossover 
region in Fig. [3jb), rather than in the clean or diffu- 
sive limits. With the typical length of CrC>2 in Q of 
L w 300 nm w ll£o> we can from Fig. [5jb) predict that 
if measurements are extended to lower temperatures, or 
a shorter junction length is used, a peak in J C (T) should 
be observed. 

In conclusion, we have presented a detailed study 
of spin-triplet supercurrents through half-metallic ferro- 
magnets. We have pointed out that ferromagnets with 
biaxial asymmetry are particularly suitable for creating 
such triplet supercurrents. We have studied in detail the 
role of disorder and found that the mechanism we propose 
is robust. We have clarified the symmetries of the pairing 
amplitudes that determine the Josephson current, and 
have made predictions that can be tested experimentally. 
In particular, there is a robust peak in the temperature 
dependence of the critical Josephson current that signals 
the involvement of odd frequency pairing amplitudes. 

METHODS 

We obtain the Josephson current as function of impu- 
rity concentration, temperature and junction length us- 
ing the quasiclassical Green's functions technique 0, 0] . 
The Green's functions g(pj?,R, e n ) depend on the spa- 
tial coordinate R, Matsubara energy e n = (2n + 1)ttT, 
and the momentum direction on the Fermi surface p F . 
In the superconductors, the propagator g is a 4x4 matrix 
in combined spin and particle-hole space, 

§S = f _9s + g t <7 (fs + ft ■ °)i>°~v I ( 7 ) 

V ^ s + ft ' a *) i(J y gs + gt-v* J ' 

where f s and f t are singlet and triplet pairing ampli- 
tudes, g s and g t are spin scalar and spin vector parts 
of the diagonal Green function, and the vector er = 
(cr x , <7 y ,o~ z ) is composed of Pauli spin matrices. The hole 
amplitudes are related the particle amplitudes by the 
symmetry f(p F ,e n ) = f(-p F ,e n )*. In the half metal, 
only conduction electrons with spin up exist, and the 
propagator is a 2x2 matrix in particle-hole space, 

g H = ( 9 ~ TT /tt ) ■ (8) 

The propagators are connected at the interfaces via the 
scattering matrices given in Eq. (fTj). 



6 



The transport equation governing the supercurrent in 
the heterostructure is given by the Eilenberger equation 
for the propagator g. Impurities are treated in the Born 
approximation using a life time of quasiparticlcs in 
the superconductor, and a life time tr in the half metal. 
The corresponding mean free paths are is — vsts, Ih — 
vhtr, with the Fermi velocities vs and vr in the two 
materials. The equation of motion for the 4x4 Green's 
function in the superconductors reads, 



iv s n d x g 



1 



2ttts 



(9), 9 



6, (9) 



where /1 = cos(O p ), 9 p is the angle between the Fermi 
velocity and the i-axis, T3 is the third Pauli matrix in 
particle-hole space, and A = Alia y is the singlet order 
parameter. The average (• • • ) = J iSSS^zjdVs . . . j s over 
all momentum directions. There is an analogous equation 
for the 2x2 Green's function in the half metal, 



iv H [i d x g 



1 



2-kth 



(9), 9 



6. (10) 



Eqs. ([9f- (fT0|) are supplemented with the normalization 
condition g 2 — — 7r 2 l. 

We linearize the above equations for small triplet com- 
ponents in the superconductor (ft ) and in the half metal 
(/ft). The normalization condition is then used to elim- 
inate the diagonal part of g in favor of a coupled set 
of equations for /-functions with positive and negative 
momentum directions. We decouple the system of differ- 
ential equations by introducing the new triplet functions 
ipo and ip in Eqs. ©-©■ 

The solutions for the functions Voj") appearing in the 
ansatz ((2]) for the superconductors, are given by 



Ss f 
V02( a; ) = -tt b 02{x) + dx 

1 Jx 2 



Xl Ki (x x') 

& ^nr^iriiV)) (11) 

00 



, K 2 (x,x') 
2Ws 



(v&OO) (12) 



with s e = sgn(e„), B 03 (x) = e -\*-*j\/esM > Kj (x,x') = 

e -\x-x'\/£s\n\ _|_ e — JCjI+Isk— a=al)/Cs|Ml) ; anc [ £ s — 

v s /(2n n +T s 1 ) withft„ = yje 2 n + |A| 2 . The momentum- 
antisymmetric parts are obtained by using the identity 
<Poj = -VSe£sd x (f s 0j . 

In the half metal the solutions for <p S j' a are 



1 - e - 2i /«ffM V 2 



(13) 



with £ H = v H /(2\e n \ + t^- 1 ), L = x 2 - Xi, Bi(x) = 
Vhfa) (e-^ 2 -^"^ + e -t L+x - Xl Vt''M), K(x,x') = 



e —\x—x'\/£n\l*\ _|_ e -(2£-|a;— 3t'|)/£h I A* I _|_ e —(x+x'—2xi)/£H\l*\ + 
e -{2x 2 -x-x )/€-»■ _ xhe momentum-antisymmetric parts 
are obtained by using the identity ip® — — /is s £#9 x <^. 

The integral equations (fTTjl - fHij) for (<p s (x)) are solved 
by replacing tp on a spatial grid with a piecewise linear 
function. Exact integration of the resulting expressions 
reduces the problem to a simple matrix inversion. The 
angular averages can be performed analytically and lead 
to exponential integrals. This procedure is necessary be- 
cause the integration kernel decays on a different length 
scale compared with ip in the diffusive limit. 

The current density in the half metal is given by the 
diagonal Green's function, 



J x (x) = ev H N H T y^(/ig-[-|-(/i, e„, a;)) 



(14) 



where Nu is the density of states at the Fermi level for 
the conducting spin band in the half metal, and e is the 
electronic charge. It can be shown from the transport 
equations © and ([IT))) that the current density J x does 
in fact not depend on the spatial coordinate, in agree- 
ment with the continuity equation that expresses particle 
conservation. 

Using the normalization condition in the half metal, 
g 2 | = — 7r 2 + /tt/tt, for small triplet amplitudes, one ob- 
tains in leading order g-\-\ = — i7rsgn(e„)(l — f^f^/2-K 2 ), 
leading to 



J x 



ievijN 



2n 



(15) 

It is instructive to decompose the anomalous propaga- 
tors into their symmetric and antisymmetric part f*^ a 
with respect to \x. Doing this and using the fundamental 
symmetries /f r (-£ n ) = -/f T (e n ), /f T (-£n) = /n( £ «) 
we arrive at 



Jx = 



2eV ^Tj2 r^MM/TVf*)- (I 6 ) 

n e n >0 J 



Substitution of Eq. © into Eq. (QUI leads to Eqs. (g))-©. 

To study the symmetry properties of the Josephson 
current, we expand the pairing amplitudes in Legen- 
dre polynomials. Writing ^ 3 (/i) = Aj{ij)tpj{ij), the ex- 
pansion is Tp~j{ii) = J2h^o ^(/•OWj./j with components 
Tpjj = {21 + l){Pi{fj,)lp(fi)}. Using that 



E 

1=0 



(-i)'q + i) 

(2Z + l)(2Z + 3) 



{fz.lfl.l+l - <Pl,l<P2,l+l) )( 17 ) 



we can bring Eq. ([5]) into the form J c = Jc-,i,i+i- 
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SUPPLEMENTARY MATERIAL 

In Supplementary Fig. [T]we show solutions of the inte- 
gral equations in the superconductor, Eqs. (fTl"j) - (fT!?]) . and 
in the half metal, Eq. (fT5|) , for several impurity concen- 
trations ranging from the ballistic limit to the diffusive 
limit. In the left panels we vary the superconducting 
mean free path is, and in the right panels the half metal 
mean free path £h for fixed is- In the inset we also show 
the induced spin polarization in the superconductors. It 
is calculated from the the diagonal part of the Green's 
function, given by 



|A| 2 ^- 



sgn(e n )cp s 0j (x). 



(18) 



There is, consequently, a surface spin polarization in 
the superconductor, that in the clean limit is given by 
(g ZJ ( x )) = -^|A| 2 (^e~ |a; ~ Xjl/?sH }/2(4 + |A| 2 )- The 
induced spin magnetization is then 



6M z (x) =2 flB N s Ty(g zj (x)) 



(19) 



where (1b is the Bohr magneton and N$ is the density of 
states in the superconductor. 
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FIG. 1: (Supplementary Figure) Spatial dependences of pair- 
ing amplitudes in the superconductor and in the half metal 
(length L = 2£o, only left half shown) for a 7r -junction. The 
odd-frequency s-wave triplet amplitude is shown in the up- 
per panel, the even-frequency p-wave triplet amplitude in the 
lower panel. In the inset we show the induced spin polariza- 
tion of quasiparticles near the interface in the superconductor. 
The triplet amplitudes in the superconductor are m = with 
respect to the superconducting quantization axis, and in the 
half metal are equal-spin m — 1 amplitudes with respect to 
the magnetization axis M. The quantization axis in the su- 
perconductor can be slightly misaligned with respect to that 
in the half metal as result of the breaking of spin-rotation 
symmetry around M at the interface. 



